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The cross section for dissociative recombination of BeH + is calculated by solution of the time- 
dependent Schrodinger equation in the local complex potential approximation. The effects of cou- 
plings between resonant states and the Rydberg states converging to the ground state of the ion 
are studied. The relevant potentials, couplings and autoionization widths are extracted using ab 
initio electron scattering and structure calculations, followed by a diabatization procedure. The 
calculated cross sections shows a sizable magnitude at low energy, followed by a high-energy peak 
centered around 1 eV. The electronic couplings between the neutral states induce oscillations in the 
cross section. Analytical forms for the cross sections at low collision energies are given. 
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I. INTRODUCTION 

The ITER fusion reactor planned for Cadarache, 
France, is the next major step in fusion research. It is 
presently scheduled to be taken into operation in 2016. A 
critical choice in the construction of ITER is the plasma 
facing material. There is not a single choice which will 
fulfil all possible criteria, and research and development 
is required in order to make an optimal choice. Much ef- 
forts have been devoted to obtain reaction cross sections 
for atoms and molecules present in the edge plasma and 
the divertor region, not least through the Coordinated 
Research Projects organised by IAEA [1]. Beryllium is 
a primary choice for the first wall because of its low ero- 
sion and atomic number. Even with a low degree of ero- 
sion, however, beryllium will enter the plasma and it is 
therefore of importance to include surface and volume 
processes for beryllium in various forms in modelling the 
edge plasma. Beryllium-containing hydrides will form in 
the edge plasma, and these molecules will be ionized by 
electrons. The question posed in this paper is how effec- 
tively the BeH + ions and its isotopologs are destroyed by 
electrons. The dominant process destroying BeH + in the 
edge plasma is dissociative recombination. This process 
has never, as far as we know, been studied experimen- 
tally, because of the toxicity of beryllium. Thus, there is 
a need for theoretical support to determine the cross sec- 
tion for this reaction. In this paper we present the first 
calculation of the dissociative recombination of BeH + . 

The relevant potential curves and autoionization 
widths are determined using ab initio electron scattering 
and structure calculations. This is followed by a diabati- 
zation procedure, where the electronic couplings between 
the neutral states can be determined. This is depicted in 
section [n] The nuclear dynamics is then studied using a 
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wave packet propagation method described in section HTTl 
The results are presented in section lTVl Atomic units are 
used throughout the paper unless otherwise stated. 



II. CALCULATION OF RELEVANT 
ELECTRONIC STATES 

The potential energy curves of the electronic states rel- 
evant for dissociative recombination of BeH + are deter- 
mined by combining structure calculations with electron 
scattering calculations. For the structure calculations, 
the Multi-Reference Configuration Interaction (MRCI) 
technique is used to determine the neutral electronically 
bound adiabatic states situated below the potential en- 
ergy curve of the ground state (A 1 S + ) of the BcH + ion. 

In a quasidiabatic representation, some of these neutral 
states will cross the ionic ground state and become res- 
onant states, i.e., they will couple to the ionization con- 
tinuum and have a finite probability for autoionization. 
In order to determine the energies and autoionization 
widths of the resonant states, the Complex Kohn Varia- 
tional method Q is used. In order to obtain the quasidi- 
abatic potentials, the scattering data is combined with 
a diabatization procedure of the electronically bound 
states. 

Finally, the calculated potential energy curves, au- 
toionization widths and electronic couplings between the 
neutral states can be used in order to carry out calcula- 
tions on the nuclear dynamics. 



A. Structure calculations 

The adiabatic potential energy curve of the BeH + 
ground state (A 1 S + ) as well as several excited states 
of BeH of 2 E + , 2 II, and 2 A symmetries are calculated 
using the MRCI technique. 



2 



For both the structure and scattering calculations, nat- 
ural orbitals are used. The natural orbitals are deter- 
mined using a MRCI calculations on the BeH + ground 
state. In this calculation, the reference space consists of 
the lcr, 2a, 3cr, Itt, 4<t orbitals and single and double exci- 
tations from the reference configurations into the virtual 
orbitals are included. The natural orbitals are calculated 
using a basis set for the H atom of (4s, lp) contracted to 
[3s, lp] 3, while for Be, a basis set of (10s,6p,ld) con- 
tracted into [3s,3p,ld] is used [H, The natural orbitals 
are then further expanded with diffuse orbitals in order 
to accurately describe the Rydberg character of some of 
the electronic states. The H-orbitals are augmented with 
(2s,2p,ld) orbitals, and the Be-orbitals with (4s,lp,ld). 

The MRCI calculations on the ionic and neutral ex- 
cited electronic states are carried out using a refer- 
ence space consisting of the la, 2a, 3<r, lit, 4a, 5a orbitals. 
These natural orbitals all have an occupation number 
greater than approximatly 0.002. 

In addition single excitations out of the reference space 
are included. The calculations are carried out in C% v 
symmetry. For the neutral system, the MRCI calcula- 
tions consist of about 5500 configurations in A\ symme- 
try, 3700 in Ai and 4600 in By. For all symmetries, the 30 
lowest roots are calculated. The potential energy curves 
are calculated for internuclear distances ranging from 1 
ao to 10 ao- 

The adiabatic potential energy curves situated below 
the ionic ground state potential are displayed in Fig. Q] 
In (a) we show the electronic states of BcH of 2 £ + sym- 
metry and in (b) we show the corresponding states of 
2 FI symmetry. The dotted curve is the potential energy 
curve of the ionic ground state. The dashed (red online) 
curves included in the figure are the very accurate poten- 
tials of BeH recently calculated by Pitarch-Ruiz et al. @ 
using a full configuration interaction calculation with a 
large basis set. As can be seen in the figure, the potential 
energy curves calculated by Pitarch-Ruiz et al. are lower 
in absolute energy. However, our potentials show a simi- 
lar dependence as a function of the internuclear distance. 
In the present calculation, our structure calculations are 
limited by the scattering calculations we perform in the 
next step (see section iHBj) . In order to describe the res- 
onant and electronically bound states at the same level 
of theory, it is important to use the same MRCI wave 
function for the target ion in both the structure and the 
scattering calculations. In Fig. [1] we do not display any 
2 A states. Our basis set is not good enough to represent 
any Rydberg 2 A states. We only obtain one resonant 
state of this symmetry. The higher excited states of 2 A 
symmetry have threshold energies above 1.05 eV rela- 
tive to the ground vibrational level of the ion and these 
states will not be important for the total cross section of 
dissociative recombination at low collision energies. 
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FIG. 1: Adiabatic potential energy curves of BeH of (a) 2 E+ 
symmetry, and (b) 2 II symmetry. The black dotted curve is 
the ground state (X 1 ^) of the BeH + ion. The black full 
curves are the potentials calculated here, while the dashed 
(red online) curves are potentials calculated by Pitarch-Ruiz 
et al. @| 



B. Electron scattering calculations 

By using the Complex Kohn variational method 0], 
the energy positions and autoionization widths of the res- 
onant states are determined. As a target wave function 
4>i(ri, ...,rN), the same MRCI wave function for the ion 
as used in the structure calculations is applied. In the 
Complex Kohn variational method, the N + 1 electron 
trial wave function is written as a 

* = ^A[$ i (r 1 ...r N )F i (r N+1 )]+^2d fl ®i l (r 1 ...r N+1 ). 

i fi 

The first sum is denoted as the P-space portion of the 
wave function and runs over the energetically open tar- 
get states. In the case of BeH, only one channel was 
open. The function Fi (rjv+i) is the one-electron wave 
function describing the scattered electron, and A is an 
anti-symmetrization operator for the electronic coordi- 
nates. The second term, denoted as the Q-space portion 
of the wave function, contains the functions which 
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are square-integrable N + 1 configuration state functions 
(CSFs) which are used to describe short-range correla- 
tions and the effects of closed channels. The advantage 
of using natural orbitals is that the orbital space used 
to generate these states is kept manageable small. The 
one-electron scattering wave function Fi is in the case of 
electron-ion scattering, further expanded as 

Fi(r) = £>}<^(r) + (2) 

i 

[fl (hr)5u S mmo +Tl lomm J+(k l r)] Yi m (r)/r. 

Here <fij(r) are a set of square-integrable functions, 
are the outgoing- and incoming Coulomb functions and 
Yi m are spherical Harmonics. Angular momenta up to 
I = 6 and \m\ — 4 are included in the calculation. 

By inserting the trial wave function into the Complex 
Kohn Functional [j| , the unknown coefficients in the trial 
wave function can be optimized. Then also the T-matrix 
(T[ lommo ) for elastic scattering is obtained and by fitting 
the eigenphase sum of the T-matrix to a Breit-Wigner 
form |7j, the energy positions and autoionization widths 
of the resonant states can be calculated. These electron- 
scattering calculations are carried out for a fixed geom- 
etry of the target ion. For BeH, the five lowest reso- 
nant states of 2 X + and 2 n symmetries as well as the low- 
est resonant state of 2 A symmetry are calculated. The 
electron scattering calculations are carried out for 1.0 ao 
< R < 3.0 ao with steps of 0.25 ao- 

C. Diabatization 

For internuclear distances smaller than 5.0 ao, the elec- 
tronic ground state of BcH + has the dominant configu- 
ration (la) 2 (2a) 2 . All Rydberg states converging to this 
ionic core have this configuration plus an extra electron 
in an outer, Rydberg-like orbital. The resonant states 
are in general Rydberg states converging to excited ionic 
cores. For these states the (2a) orbital is typically singly 
occupied. It is therefore relative easy to use the config- 
uration interaction coefficients of the MRCI wave func- 
tions to follow the resonant states when they cross the 
Rydberg manifold situated below the ionic ground state 
potential. This is done in order to obtain an "estimate" 
of the quasidiabatic potential energy curves. 

It should be mentioned that for R > 5.0 ao, the 
ionic ground state changes dominant configuration to 
(l<r) 2 (2cr) 1 (3(T) 1 . For these large internuclear distances, 
the Rydberg states will then have similar configurations 
as the resonant states and the quasidiabatic potentials 
are simply estimated by the shape of the adiabatic po- 
tentials. 

This technique, however, will only provide us with di- 
agonal element of the diabatic potential matrix and not 
the electronic couplings between the diabatic states. In 
order to obtain these couplings a diabatization proce- 
dure is applied. This diabatization procedure was very 



recently developed by us for determining the quasidia- 
batic potential matrix for the HF molecule @. There 
is a unitary transformation matrix T that will trans- 
form the adiabatic potential matrix [with matrix ele- 
ments Vfi — V i a (R)Sij] into the diabatic potential matrix 

V d = TV a T-\ (3) 

In this diabatization procedure, we assume that the 
transformation matrix can be written as a product of 
matrices that describe successive 2x2 rotations among 
the adiabatic states T = . . . T2T1. The matrices are 
of the form (in the case of a rotation among state 1 and 
2) 

/ cos(7i) - sin(7i) . . . \ 
sin(7 l ) cos(7i) 0... 

T 10. (4) 

1 ... K ' 

V : ! ■■) 

where we assume that the rotational angles have the an- 
alytical form 

7i = ^[l + tanh(a i ( J R- J R i ))]. (5) 

We then perform an optimization procedure, where the 
unknown parameters and Ri of the rotational angles 
are determined in order to optimize the agreement be- 
tween the diagonal elements of the diabatic potential ma- 
trix in Eq. ([3]) and the estimated quasidiabatic potentials 
obtained using the Cl-coefficients as described above. 
When the transformation matrix is optimized, not only 
the diagonal diabatic potential curves are obtained, but 
also the electronic couplings between the neutral states. 
In Fig. [5] we show the diabatic potential curves of BeH 
obtained using this procedure. 

The resonant states situated above the ionic ground 
state are calculated using the electron scattering calcu- 
lations. These resonant states will then cross some of 
the Rydberg states situated below the ionic potential. 
Where these curves cross, they will couple by electronic 
couplings. In Fig. [31 we show examples of the electronic 
coupling elements between the 2 S + neutral states marked 
with the thicker curves in Fig. [2] with the other states 
of the same symmetry. We label the diabatic electronic 
states from 1 to n in energy order from lowest to highest 
at the internuclear distance 1.0 ao. It should be noted 
that the electronic couplings are localized to the regions 
of avoided crossings. The couplings between the electron- 
ically bound state 2 and the higher states [see Fig [3] (a)] 
have the double peak structure since the compact state 
2 crosses several of the Rydberg states twice. This state 
does not, however, cross the ionic ground state. 

Similar couplings, but now for the 2 n manifold are 
displayed in Fig 2] 



4 




(a) 0.02 



a 0.01 - 



8 0.005 





(C)-14.6 



R(a n ) 




R(a n ) 



FIG. 2: Quasidiabatic potentials of BeH of (a) 2 S + , (b) 2 II 
and (c) 2 A symmetries. The solid (red online) curves are 
the resonant states while the electronically bound states are 
dashed (black) curves. Also the electronic ground state of the 
BeH + is displayed with the dotted curve. 



III. NUCLEAR DYNAMICS 



The cross section for dissociative recombination of 
BeH + is calculated by propagating wave packets on the 
neutral states. The electron capture will induce the ini- 
tial condition for the wave packets on the resonant state 
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FIG. 3: Some of the electronic couplings between diabatic 



states of 2 E + symmetry. 



*<(t = 0,iZ) 



2tt 



Xv(R)- 



(6) 



Here Ti is the autoionzation width for the resonant state 
and Xv is the vibrational wave function of the ion. The 
rotational motion of the molecule is neglected. 

The nuclear dynamics is studied by numerically inte- 
grating the time-dependent Schrodinger equation 



•V d *(t,iJ) (7) 



using a Cranck- Nicholson propagator 10]. In the present 
study the nuclear dynamics is studied using both cou- 
pled and uncoupled electronic states. For the uncoupled 
states, Vd is a diagonal matrix with the quasidiabatic po- 
tentials of the resonant states on the diagonal. Autoion- 
ization is included by letting the resonant states become 
complex when the resonant state potential is situated 
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The wave packets are propagated out into the asymp- 
totic region, where the cross section for dissociative re- 
combination is calculated [UdH] by projecting the asymp- 
totic wave packets onto the energy-normalized eigen- 
statcs of the fragments 

*i(E) = ^g\($%(R)\Mtoo,R))\ 2 . (9) 

Here g is the ratio of multiplicity between the the res- 
onant state and ionization continuum. The total cross 
section for dissociative recombination is obtained by sym- 
marizing all partial cross sections a(E) — J2i 

In the present calculations, the wave packets are prop- 
agated on a grid ranging from 0.5 ao to 300 ao with 
dR = 0.01 ao- The wave packets are propagated with 
time-steps of dt = 0.1 a.u. For the uncoupled potentials, 
the wave packets are propagated until too = 1000 a.u., 
while too — 4000 a.u. is needed for the coupled system. 



IV. RESULTS AND DISCUSSION 

As mentioned above two models are applied in the 
study of dissociative recombination of BeH + . In the 
first model, only the dynamics on the uncoupled reso- 
nant states is considered, while in the second model the 
propagation of the wave packets on the coupled poten- 
tials is explored. In both models, autoionization from 
the resonant states is included using local complex po- 
tentials. 



FIG. 4: Some of the electronic couplings between diabatic 
states of 2 II symmetry. 

above the ionic ground state 

Vf{R) = Vf{R)-i^-. (8) 

This is the so-called local approximation for treating au- 
toionization @|. The validity of this approximation is 
examined using a non-local expression for autoioniza- 
tion [ll[ and we found that even for the two lowest reso- 
nant states of BeH of 2 II symmetry, the local approxima- 
tion is valid. When the electronic couplings are included, 
the matrix in eq. ||7J) now contains both the resonant 
and electronic bound state potentials on the diagonal as 
well as the off-diagonal coupling elements. The electron 
capture directly into the Rydberg states as well as au- 
toionization out of the Rydberg states are in the present 
study neglected. The neglect of the so-called indirect 
mechnism is valid when there is a strong direct mech- 
nism [ijl, 0, HH, i.e., when there are resonant states 
with potentials crossing the ion close to its minimum. 
The indirect mechnism will only affect the cross section 
at very low collision energies and create sharp structures 
in the cross section. 



A. Uncoupled states 

As described above, five resonant states of 2 £ + sym- 
metry are included in the model. Only the two lowest 
resonant states (here labelled with and E/gg) are asso- 
ciated with asymptotic limits that arc energetically open 
at zero collision energy. The following 2 £ + states open 
up for dissociation at energies of 0.56 eV (U\q\q), 0.87 eV 
(C^ini) and 1.59 eV (C/1212) relative to the ground vibra- 
tional level of BeH + . In the following the cross section 
for electron recombination with BeH + in its ground vi- 
brational state is calculated, that is, the vibrational wave 
function of the v = level of BcH + is used in Eq. ©. 
In Fig. [5] (a), we display the cross section for dissocia- 
tive recombination, obtained when propagating the wave 
packets on the uncoupled 2 £ + states. As can be seen 
in the figure, the cross sections originating from all 2 S + 
states show a smooth onset. This reflects the fact that 
the resonant states for this symmetry do not cross the 
ionic ground state potential close to the equilibrium sep- 
aration of the ion. The energy-dependence at threshold is 
therefore determined by the capture probability which is 
related to the Franck-Condon factors between the initial 
vibrational wave function of the ion and the continuum 
wave function of the resonant states. 
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The remaining states of 2 II symmetry have threshold en- 
ergies of about 0.83 eV (for V77 and V$a) and 1.81 eV for 
Vgg. As can be seen in Fig.[5](b), where the cross section 
from the 2 I1 states are plotted using a linear-log scale, it 
is the low energy tail of the cross section from the second 
resonant state (V55) that dominates the low-energy cross 
section for this symmetry. 

Finally, in Fig. [5] (c), the cross section from the 2 A 
resonant state is displayed. This state opens up for dis- 
sociation at 0.63 eV and the cross section shows a sharp 
threshold at this energy. 

In Fig. the total cross section using a log-log scale 
calculated with uncoupled resonant states is shown. In 
the figure, also the contribution to the cross section from 
states of the different symmetries are displayed. As can 
be seen, at low energy the cross section is dominated by 
resonant states of 2 II symmetry. The sharp threshold 
around 0.63 eV, asising from the 2 A resonant state can 
be seen in the total cross section. 
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FIG. 6: Total cross section for dissociative recombination of 
BeH + , calculated using uncoupled resonant states. 



FIG. 5: Cross section for dissociative recombination of BeH , 
calculated by propagating wave packets on uncoupled reso- 
nant states of (a) 2 E + symmetry, (b) 2 II symmetry (using a 
linear-log scale) and (c) 2 A symmetry. The thin curves are 
the cross sections originating from each resonant state, while 
the thick curve is the total cross section from the states of a 
given symmetry. 



The situation is different for the resonant states of 2 II 
symmetry, where two states (labelled with V33 and V55) 
cross the ion potential at smaller distances than the equi- 
librium distance. The first resonant state (V33) has a 
threshold of about 0.63 eV, while the second state (V55) 
is energetically open for dissociation at zero collision en- 
ergy. The cross section from this state shows the typical 
X/E dependence [l6| of the cross section at low energies. 



B. Coupled states 

As mentioned above, several resonant states of BeH 
are not energetically open for dissociation at zero col- 
lision energy. We therefore wish to examine a model 
where the electronic couplings between the resonant and 
electronically bound states are included. We have pre- 
viously seen for other systems such as HF + 0, |T3| that 
the inclusion of electronic couplings between the neutral 
states may open up pathways to dissociation. Flux that 
in the uncoupled model is trapped in bound potentials 
will transfer to lower states that are open for dissocia- 
tion. 

As described above in section III CI the neutral states 
are diabatized and the electronic couplings are deter- 
mined. For the states of 2 E + symmetry, seven states 
are diabatized, i.e., the lowest resonant state and six 
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electronically bound states. In a quasidiabatic represen- 
tation, there is one compact electronically bound state 
that will cross several of the Rydberg states situated be- 
low the ionic ground state. However, this state will never 
cross the ion potential and become a resonant state. The 
ground state of BeH, X 2 £ + , is assumed not to couple to 
the higher electronic states. Due to numerical limitations 
in the optimization procedure, the electronic couplings to 
the higher resonant states of this symmetry are not cal- 
culated. The wave packets are then propagated using the 
full diabatic potential matrix. In Fig. [71 the cross section, 
calculated using the coupled 2 £+ states is shown. In the 
figure, also the cross section calculated using uncoupled 
states are dispalyed. 
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FIG. 7: Cross section calculated using the coupled states of 
BeH of 2 £ + symmetry. The thin solid curves are the cross 
sections from the resonant states, while the dashed curves are 
cross sections from the electronically bound states. The total 
cross section originating from the coupled 2 £ + state is shown 
with the thick black curve. Also the cross section calculated 
using uncoupled 2 £ + states is displayed with a thick dashed 
black curve. 

The inclusion of the electronic couplings will not dra- 
matically change the magnitude of the cross section. 
Sharp structure is created both in the partial cross sec- 
tions as well as the total cross section. This structure 
can be interpreted as Feshbach resonances created when 
the states open for dissociation are coupled to the nu- 
clear bound states. Structures in the cross section is then 
formed when the wave packets can couple from a reso- 
nant state and form a vibrationally bound state in an- 
other potential. Similar structures have been seen in the 
cross section of HF + 8] . Above the dissociation energy 
of the ion (around 2.66 eV) there are no bound vibra- 
tional levels of the Rydberg states and this explains the 
lack of structures at higher energies. Furthermore, the 
dynamics on the highest 2 £ + resonant states is studied 
using uncoupled potentials. 

When the electronic states of 2 n symmetry are dia- 
batized, eight states are included, i.e., all five resonant 
states and three electronically bound states. The low- 



est A 2 H state is not included in the process. Again, the 
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FIG. 8: Cross section calculated using the coupled states of 
BeH of 2 n symmetry. Also the cross section calculated using 
uncoupled resonant states is displayed. 

electronic couplings induce sharp oscillations in the cross 
sections, but they will not significantly effect the magni- 
tude of the cross section. 

The fact that the inclusion of the couplings do not in- 
crease the magnitude of the cross section at low energies 
may be understood by smooth onset of most of the cal- 
culated partial cross sections. Hence, it is the electron 
capture that limits the energy components of the wave 
packets captured into the resonant states. Simply, there 
is no low-energy components initiated in the resonant 
states that can find a way to dissociation by inclusion of 
electronic couplings. 

Fig. [5] shows the total cross section for dissociative 
recombination of BcH+ calculated using both the un- 
coupled and coupled models. Both models produce a 
cross section with a similar magnitude, but oscillations 
are created when the electronic couplings are included. 
The cross section has a low-energy 1/E dependence and 
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FIG. 9: Total cross section for dissociative recombination of 
BeH + , calculated using coupled and uncoupled potentials. 
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a mag nitude of about 3 x KT 14 cm 2 at E = 0.01 eV. 
This is a typical magnitude for a dissociative recombina- 
tion cross section of a diatomic molecule [HI- Further- 
more, the cross section shows a high-energy peak centered 
around 1 eV. Similar high-energy peaks have been seen in 
the cross section for dissociative recombination for sev- 
eral ionic species [llj]. The sharp threshold around 0.63 
eV in the total cross section comes from the 2 A resonant 
state. 



C. Effects of vibration 

The cross section for dissociative recombination with 
BeH + in the v = 1 vibrational level is calculated by using 
the Xv=i vibrational wave function in the initial condition 
of the wave packets [eq. (|6[)]. Furthermore, the energy- 
scale is shifted to account for the shift in energy of the 
target ion. The wave packets are then propagated on the 
coupled potentials. In Fig.[ini we compare the total cross 
section for dissociative recombination of BcH + in v = 
and v = 1. It can be seen that the shift in the energy 
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the mass has to be changed. Again, the coupled po- 
tentials are used for calculating the cross section. The 
resulting cross section is displayed in Fig. [IT] and com- 
pared with the corresponding cross section for BeH + . 
The cross section of the two isotopologues are similar 
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FIG. 11: Total cross section for dissociative recombination of 
BeD + and BeH + in their lowest v = vibrational states. 

in shape and magnitude. They both show oscilliations 
induced by the electronic couplings as well as the high- 
energy peak around 1 eV. The magnitude of the high- 
energy peak is smaller for the heavier isotopologue and 
similar behaviours have been seen for the high-energy 
peak oberserved in for example dissociative recombina- 
tion with Hg" [l9(. The cross section for electron recom- 
bination with BcD + in the vibrationally excited v = 1 
state is calculated and similar to BeH + , the positions of 
the structures in the cross section as well as the high en- 
ergy peak are shifted in the energy and the magnitude of 
the low-energy cross section is increased. 



FIG. 10: Total cross section for dissociative recombination of 
BeH + in the v = 1 and v = 1 vibrational states. 

scale causes shifts of the resonant structures as well as the 
high energy peak towards smaller energies. Furthermore, 
the cross section at low collision energies is increased and 
this can be understood by larger capture probability into 
the low lying resonant state for the vibrationally excited 
ion. 



D. Effects from isotopic substitution 

The cross section for dissociative recombination with 
BeD + in the v = vibrational state is calculated by us- 
ing the same set of potentials, widths and couplings, but 
changing the reduced mass to the mass of BeD. Both 
when the vibrational wave function of the target ion 
is calculated as well as in the wave packet propagation 



E. Analytical forms for the cross sections 

To simplify the use of the calculated cross sections in 
modelling of the fusion plasmas, the cross sections have 
been fitted to analytical forms. The cross sections calcu- 
lated using uncoupled potentials have been used for the 
fitting. These cross sections show a smooth energy de- 
pendence at collision energies below about 0.3 eV that 
easily can be fitted to functions of the form 

°m = § an) 

In table HI the values of the parameters <7o and b are given 
for the fits of the calculated cross sections for dissociative 
recombination of BeH + and BeD + in the v = and v = 1 
vibrational levels. 
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TABLE I: Parameters for the fitted form of the dissociative 
recombination cross sections at low collision energies. 



Ion 



vibrational level cr (10 16 cm 2 xeV) 



BeH+ 
BeH+ 
BeD+ 
BcD+ 



v = 
v = 1 
v = 
v = 1 



3.10 
6.90 
3.32 
4.72 



0.99 
0.93 
1.02 
1.05 



CONCLUSIONS 



The reaction is studied assuming both uncoupled and 
coupled potentials. The inclusion of the couplings cause 
oscillations in the cross section and will not significantly 
affect the magnitude of the cross section. The cross sec- 
tion shows a low-energy 1/E dependence followed by a 
high-energy peak centered around 1 eV. The cross sec- 
tion for dissociative recombination with vibrationally ex- 
cited v = 1 ions is larger at low collision energies than 
the corresponding cross section for vibrationally relaxed 
ions. The cross section for the heavier isotopologue BeD + 
is similar to the cross section for BeH + , except for the 
high-energy peak that becomes smaller in magnitude. 



For the first time, the cross section for dissociative re- 
combination of BeH + has been calculated. This cross sec- 
tion is important in the modeling of the fusion plasmas. 
However, the toxicity of BeH + makes measurements on 
the reaction extremely difficult. In the present study, 
structure and electron scattering calculations are com- 
bined to obtain qausidiabatic potentials of BeH. Using an 
optimization procedure, the electronic couplings between 
the resonant states are estimated. The dynamics of the 
reaction is then explored using a wave packet technique. 



Acknowledgments 

We are thankful for data of calculated potential curves 
send to us by Jose Sanchez-Marin. A. L. acknowledges 
support from The Swedish Research Council and A. E. O. 
acknowledges support from the National Science Foun- 
dation, Grant No. PHY-05-55401. M. L. participated 
in the IAEA Coordinated Resreach Project "Atomic and 
Molecular Data for Plasma Modelling" . 



[1] R.E.H. Clark (ed.) 2006, Summary Report of Final IAEA 
Research Co-ordination Meeting. Data For Molecular 
Processes in Edge Plasmas, INDC International Nulcear 
Data Committee, INDC(NDS)-0491, IAEA, Vienna. 

[2] T. N. Rescigno, C. W. McCurdy, A. E. Orel and B. H. 
Lengsfield III, "The Complex Kohn Variational Method 
in Computational Methods for Electron- Molecule Scatter- 
ing" eds. W. H. Huo and F. A. Gianturco (Plenum, New 
York (1995)). 

[3] T. H. Dunning Jr., J. Chem. Phys., 53 2823 (1970). 

[4] T. H. Dunning Jr. and P.J. Hay, in Methods of Electronic 
Structure Theory, vol 2 (H.F. Schaefer III, ed., Plenum 
press (1977)). 

[5] D. Feller, J. Comp. Chem., 17, 1571 (1996). 

[6] J. Pitarch-Ruiz, J. Sanchez-Marin, A. M. Velasco, and I. 
Martin, J. Chem. Phys., 129, 054310 (2008). 

[7] S. Geltman, "Topics in Atomic Collision Theory" (Aca- 
demic Press, 1997) p. 31. 

[8] J. B. Roos, A. Larson and A. E. Orel, submitted to Phys. 
Rev. A (2009). 

[9] C. W. McCurdy and J. L Turner, J. Chem. Phys., 78, 



6773 (1983). 

[10] A. Goldberg, H. M. Schey and J. L. Schwartz, Am. J. 

Phys., 35, 177 (1967). 
[11] A. Giusti-Suzor, J. N. Bardsley, and C. Derkits, Phys. 

Rev. A, 28, 682 (1983). 
[12] D. R. Bates. J. Phys. B: At. Mol. Opt. Phys., 24, 695 

(1991). 

[13] A. I. Florescu, V. Ngassam, I. F. Schneider, and A. Suzor- 
Weiner, J. Phys. B: At. Mol. Opt. Phys., 36, 1205 (2003). 

[14] V. Ngassam and A. E. Orel, Phys. Rev. A, 73, 032720 
(2006). 

[15] K. C. Kulander and E. J. Heller, J. Chem. Phys., 69, 

2439 (1978). 
[16] E. P. Wigner, Phys. Rev., 73, 1002 (1948). 
[17] J. B. Roos, A. Larson and A. E. Orel, Phys. Rev. A, 78, 

022508 (2008). 

[18] A. Al-Khalili et al, Phys. Rev. A, 68, 042702 (2003). 
[19] M. Larsson, H. Danared, A. Larson, A. Le Padellec, J. 

R. Peterson, S. Rosen, J. Semaniak, and C. Stromholm, 

Phys. Rev. Lett., 79, 395 (1997). 



